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The failure diagnosis for discrete event systems (DESs) has been given considerable attention in recent 
years. Both on-line and off-line diagnostics in the framework of DESs was first considered by Lin Feng in 
1994, and particularly an algorithm for diagnosability of DESs was presented. Motivated by some existing 
problems to be overcome in previous work, in this paper, we investigate the minimal cost algorithm for 
diagnosability of DESs. More specifically: (i) we give a generic method for judging a system's off-line 
;> diagnosability, and the complexity of this algorithm is polynomial-time; (ii) and in particular, we present an 

_ algorithm of how to search for the minimal set in all observable event sets, whereas the previous algorithm 
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may find non-minimal one. 
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I. Introduction 

WITH man-made systems becoming more and more complex, detecting and locating component 
failure is not a simple task. Therefore, there is a strong need for a systematic study of diagnostic 
problems and diagnosability issues [30]. As an important kind of man-made systems, discrete event system 
(DES) is a dynamical system whose state space is discrete and whose states can only change as a result of 
asynchronous occurrence of instantaneous events over time [2]. Up to now, DESs have been successfully 
applied to provide a formal treatment of many technological and engineering systems [3, 5, 16]. Naturally, 
the diagnosability of DESs is of theoretical and practical importance. 

Actually, diagnosability of DESs has received extensive attention in recent years ( for example, [1,4, 
6-9, 11-15, 17-32]). Especially, in [15], the definitions of "off-line" diagnosability and "on-line" diagnos- 
ability were introduced, and both "off-line" diagnostic algorithm and "on-line" diagnostic algorithm were 
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significantly established in the framework of DESs. However, the algorithms presented in [15] have some 
shortcomings: 1) the computational complexity of "off-line" diagnostic algorithm is exponential in general; 
2) and "off-line" diagnostic algorithm could not find the minimal one in observable events sets (OESs), 
and the algorithm of how to inspect an automaton being diagnosable was not yet given. Motivated by these 
issues, our goal in this paper is to solve these problems. 

The remainder of the paper is organized as follows. In Section JIJ we first introduce a general framework 
of diagnosability of DESs, and then explain the lost of Algorithm 1 in [15]. In Section JIIJ the definition 
of "off-line" diagnostics is first provided, and we then present a polynomial-time algorithm to realize it. In 
Section [TV] we demonstrate the principle of finding the minimal set in an OES, and particularly, present 
our new algorithm to realize it; Section [V] provides two examples to illustrate these algorithms in Sections 
irTTI and |IV] Finally some remarks are made in Section VI to conclude the paper. 

II. Preliminaries 
A. A general framework for automata and diagnostics 

1 ) DFAs: A deterministic finite automaton (DFA) can be formally defined as a 5-tuple (Q, E, 5, qo, F), 
where Q is a finite set of states, E is the input alphabet, 5 : Q x E — > Q is the transition function, qo € Q is 
the starting state, and F C Q is a set of accepting states. Operation of the DFA begins at qo, an d movement 
from state to state is governed by the transition function 5. S must be defined for every possible state in Q 
and every possible symbol in E. 

A DFA can be represented visually as a directed graph. Circular vertices denote states, and the set of 
directed edges, labeled by symbols in E, denotes 5. The transition function takes the first symbol of the 
input string, and after the transition this first symbol is removed. If the input string is e (the empty string), 
then the operation of the DFA is halted. If the final state when the DFA halts is in F, then the DFA can 
be said to have accepted the input string it was originally given. The starting state go is usually denoted by 
an arrow pointing to it that points from no other vertex. States in F are usually denoted by double circles. 

DFAs recognizes regular languages, and can be used to test whether any string is in the language it 
recognizes. As it is known, DFAs have been used to model DESs [2]. In the following, we use DFA to 
represent a DES. 

2) Diagnosability of Discrete Event Systems: We model the system to be diagnosed as a pair G = (M, E c ). 
The first component M denotes a nondeterministic Mealy automaton: 

M = (E, Q, Y, 6, h) 
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where E is the set of finite events; Q is the set of finite states; Y is the output alphabet space; 5 : E x Q — > l9 
is the state transition function. 5(a, q) gives the set of possible next states if a occurs at q; and h : T,xQ —>■ Y 
is the output function, h(a, q) is the observed output when a occurs at q. The second component E c C E is 
the set of controllable events, where the controllability of events is interpreted in a strong sense: a controllable 
event can be made to occur if physically possible. 

States of the system describe the conditions of its components. Therefore, to diagnose a failure is to identify 
which state or set of states the system belongs to. Thus, depending on the requirements on diagnostics, we 
partition the state space Q into disjoint subset (cells) and denote the desired partition by T. The state in the 
same cell are viewed as equivalent as far as failures under consideration are concerned. The model is rather 
general since we do not put any restrictions on T. 

3) Some notations: For convenience, we give some notations. Let 

M = {0-1,02, • • • ,o-„}, 
where o"j, i = 1, . . . , n, are the observed events, and the cost of M is denoted by 

C(M) = {c(a 1 ),c(a 2 ),...,c(a n )}. 

where c{ai) means the cost of observes event ai,i = 1, ... ,71, respectively. Without loss of generality, we 
suppose that 

c(o-i) > c(o- 2 ) > > c(a n ). 
Let E € OES. We denote C7(E ) = c(o"«), which presents the cost of E G , and 

minL(Yi ) = min{i : Oi € E G }. 

An important problem is how to find the smallest observable event set that makes G diagnosable for a 
given partition T. In order to solve this problem, we define the set of all observable event sets (OESs) that 
ensure the diagnosability of the system as: 

OES{T) = {E C E: G is diagnosable with respect to E Q and T}. 

B. Lost of Algorithm 1 of [15] 

1) Algorithm 1 of [15]: In order to remove events one by one in the given order until the diagnosability 
of the system is no longer ensured, Algorithm 12.11 (Fig. [IT) was presented in [15]. 

However, Algorithm 12.11 has some shortcomings, we will illustrated them in next subsection. 
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Algorithm 2.1: (minOES) 

Input: Read G = (M, E c ), M = (E, Q, Y, 6, h),T; 
\mtisAwatiow.minOES := E; 
Removal: For i = 1 to n do 
begin minOES := minOES \ {<?i}; 
if G is not diagnosable with respect to 
minOES an T then 

minOES := mmO-BiSUjcrj}; 
end; 

Output:Return minOES; 

Fig. 1. Algorithm 1 of [15]. 

2) Losf Example of Algorithm 1 of [15]: In fact, the "off-line" diagnostic algorithm (Algorithm 12. II ) could 
not find the minimal one in observable events sets. For example, let 

M = {0-1,0-2,0-3,0-4,0-5}, 

and the cost of M is 

C(M) = {13,9,7,5,2}. 

OES(T) = {£i}U{£ 2 }U{E CM:£iCS, or, E 2 C X}, where ^ = {ct 2 , cj 5 }, and E 2 = {^3, ^4, cr 5 }- 
G is diagnosable with respect to given T and an element of OES(T). 

If we use Algorithm 12. II we can get the minOES = E2, the cost of E2 is 7 + 5 + 2 = 14. However, the 
cost of £1 is 9 + 2 = 11, which is less than the cost of E 2 . Therefore, the minOES is not the minimal 
cost of OES(T). 

III. Off-line Diagnostics 

Off-line diagnostics means that diagnosis is performed when the system is not in normal operation [15]. 
For example, what a mechanic does to an automobile in a repair shop can be viewed as off-line diagnostics. 
In order to perform off-line diagnostics, one can "open" the system, access the inside, do various tests, and 
measure responses that may not be available from the system outputs. In fact, during off-line diagnostics, 
the system is not actually in operation. Therefore, the failure status of system components will not change, 
unless such changes are made in purpose. So tests can be designed with great flexibility and the order of 
testing is not critical as far as diagnosability is concerned. 



h(a, q) 
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A. Off-line Diagnostics [15] 

For off-line diagnostics, we specialize the model introduced in the previous section by assuming that the 
output are events observed. That is, Y = S , where S Q C £ is the set of observable events and the output 
map h : £ x Q — > S is a projection defined as: 

a if <7 € S D , 
e otherwise, 
where e is the empty string. 

As it was discussed before, in off-line diagnostics all events are assumed to be controllable. Therefore, 
S c = E. Since the failure states of system components will not change, information derived from all the 
test outputs are updated and relevant. 

During off-line diagnostics, if an event a G S c is observed, then the possible state of the system is: 

Q(a) ={q£Q: (3q' G Q)S(a, q') = q}. (1) 

Hence, we know every state of the system is in either Q(a) or Q — Q(a) after observing a. That is, each 
observable event partitions the state space into: 

T a = {Q(a),Q-Q{a)}. (2) 

Since there is not restriction on the tests performed in off-line diagnostics, we can observe all observable 
events that are physically possible and then determine which states the system is in. If this information is 
sufficient for us to determine which component is broken (i.e., which cell of T the system is in), then we 
say the system is off-line diagnosable. Formally: 

Definition 3.1: G is said to be off-line diagnosable with respect to T if 

f\T a <T (3) 

<TgE 

where A denotes conjunction and < means "is finer than". 

Clearly, diagnosability depends on both the observable event set S and the desired partition T. 

B. An algorithm for off-line diagnosability 

In Section [Till we have introduced the definition of "off-line" diagnostics (see Definition 13. H and equation 
© in Section Hill ). In equation ©, the right part T is given by the system. Now we must first calculate 
the left part /\ T a , where T a is given in equation © of Section Hill and Q(a) is given in equation (fl} of 
Section |III] From these two equations, given an element a, for every element of Q, it must be in Q(a) or 
not in Q{c?) (in Q \ Q(a)). So we can use one bit to identify every element of Q in Q(a) or not in Q(a) 
(i.e, 1 for elements in Q(a) and for elements not in Q(a)). Now we give algorithms to realize them. 
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Algorithm 3.2: (QC). 
Input: 5,a,Q,q; 

Initialization: Set m = \Q\,QC := 
False; 

Judge: for i = 1 to m do 

if S(a, qi) == q then 
QC := True,break; 
Output: Return QC; 

Fig. 2. Algorithm: whether a state q in Q(u) or not, that is S(a, Q) = q or <$(cr, Q) ^ q. 

Algorithm 3.3: (TQC). 
Input: 5, S , Q; 

Initialization: Set m — \Q\,n — |S |,Sj = 
0(j = l..m) 

Intersection: for i = 1 to ?i do 

for j = 1 to m do 

if <5(<Tj, Q) = ^(Algorithm U!2b then 
Sj | = (l«(i-1)); 
Output: Return Sj(j — l..m); 
Fig. 3. Algorithm: /\ T a . 

7,) Algorithm: Algorithm 13.21 (Fig. [2]) gives whether a state q in Q(a) or not, that is 5(a,Q) = q or 
6(<r,Q)?q. 

Algorithm 13.31 (Fig. [3l gives the calculation of /\ T a . 

o-eE 

In Algorithm 13.31 all the elements of F (F G /\ T^) have the same value Sj, since they have the same 

crGS 

operation in Algorithm 13.31 And /\ T a < T means that every element of /\ T a is the subset of G 

o-es ' o-es 

(G 6 T). The reverse proposition means that there exists an element of f\ T a , not all of its elements are 

the elements of G(G G T). From this, we have Algorithm 13.41 (Fig. HJ). 

2) Algorithm Complexity: In Algorithm 13.21 in "Judge" recycle, the bad time is m. So the time complexity 
of Algorithm 13.21 is 0(m). 

In Algorithm 13.31 in "Intersection" recycle, it has two loops, the complexity of first line is 0(n); the 
complexity of second line is 0(m). In third line, it calls the Algorithm 13.21 so the bad time is 0(m); and 
then the total complexity in "Intersection" recycle is 0(m 2 n). Therefore, the time complexity of Algorithm 
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Algorithm 3.4: (OFD). 






Input: S, E , Q, T; 






Initialization: Set OFD := 


- True; 


Diagnosing: Get sj(j = 


l..m 


) from 


Algorithm H31 






Applied Quicksort Algorithm to 




l..m); 






for j = 1 to m — 1 do 






if (sj == sj + i) then 






begin Find Ti g 


T, 




S.t.Oj £ Ti, 






if CTj+iGTi then 






OFL> := False; 






end; 






Output: Return OFD; 







Fig. 4. Algorithm: /\ T a <T. 



[331 is 0(m 2 n). 

In Algorithm 13.41 in "Diagnosing" recycle, it first calls the Algorithm 13.31 the time complexity is 
0(m 2 n); then it calls the Quicksort Algorithm, the bad time complexity is 0(m 2 ); for the other lines, it has 
one loop, the total complexity is 0(m 2 ). In conclusion, the time complexity of Algorithm 13.41 is 0{m 2 n). 

IV. New Algorithm for finding the minimal one in OESs 
A. Finding the minimal one in OESs 

We would like to find a minimal element in OES(T) as follows. 

Proposition 4.1: If OES(T) is not null, then the minimal elements of OES(T) exist, but may not be 
unique. 

Proof: The proof of the existence of minimal elements is straightforward. Since S is finite, 2 s is a 
finite set. Notice that S Q C S, therefore, S Q is an element of 2 s , and the elements of OES(T) are finite. 
As a result, there exists a minimal element in OES(T) . 
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The following example shows that the minimal elements of OES{T) may not be unique. Let 

£ = {a,/3,7} 

Q = {qi,Q2} 

S(a,qi) = {q 2 } 
S(P,<b) = til} 
<*(7>9i) = {92} 
5(<j,q) = otherwise, 

and 

T = {{ qi },{q 2 }}. 

Obviously, {a}, {(3} and {7} are minimal elements of OES{T). ■ 
From Proposition 14.11 in Section Hill we conclude that we may be able to find more than one set of 
observable events, and each set is minimal in the sense that removing any event from the set will make 
the system not diagnosable. Practically, we can find a cost-effective minimal observable event set by first 
ordering the events in terms of the difficulty (and hence cost) in detection. This directly gives the Algorithml 
of [15](Fig. H] Algorithm [27Tb . 

B. New Algorithm 

From the Example in Section III-B.2I we know that the minOES given by Algorithm 12.11 is not the 
minimal cost one. Therefore, we will modify Algorithm 12.11 to find the minimal cost one in this subsection. 

Proposition 4.2: By Algorithm 12. II we get the minOES, whose cost is C(minOES) and whose minimal 
label is minL(minOES). If there exists an S D G OES(T), with C(E C ) < C(minOES), then we have 

minL(Y, ) < minL(minOES). 
Proof: If mmI(E ) > minL(minOES). Set L = minL(minOES) is the minimal index of set 
minOES. In Algorithm 12.11 when / = L, G is diagnosable with respect to minOES and T, and the next 
step of Algorithm 12.11 is not executed. So ai^minOES, and then L ^ minL(minO E S) . Consequently, 
minL(Y^ ) < minL{minOES). ■ 
Now we present a new algorithm 14.31 (Fig. [5]) to find the minimal cost one. 

C. Necessary Element 

Definition 4.4: (Necessary Element) Suppose S G OES(T). If cr; e S, but E \ {ai}eOES{T), then we 
call o~i necessary element with respect to T. 
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Algorithm 4.3: (MMOES). 

Input: G = (M,E C ), M 
(E,Q,Y,S,h),T, the order £ 
{ax, o"2, . . . , er„}; the cost of S, C(E) = 
{c(cri),c(cr 2 ), ■ ■ . ,c(ct„)}; 

Initialization: Get minOES by 

Algorithm IZTl 

Set (mS = minL(minOES), 

cmS = C(minOES); 

Set i? = {S C £ : minL(Y. ) > 

lmS,C(Y, ) < cmS}; 

Set ng = \H\, and i? = 

{Hi,H 2 , ■ ■ ■ , H ng }; 

Testing diagnosability: for i = 1 to ng 

do 

begin if (G is diagnosable with 
respect to B l an T) AND 
(C(Hi) < cmS) then 

minOES = 

Hi,cmS = C{Hi); 
end; 

Output: Return minOES; 
Fig. 5. AlgorithmModify MinOES. 

Proposition 4.5: If S e OES(T), and S„ C F, then F G OES(T). 
Proof: Because 

A T - < A T - ^ T ' 

the proposition holds true. 

Proposition 4.6: If Oi is a necessary element, then for any S D € OES(T), crj G S . 

Proof: (proof by contradiction) If the theorem is not true, there exists an S G € OES(T), with (TjGS 
Therefore, S G C S \ {crj. And <7j is a necessary element, £ \ {<Ti}G0^5(T). So S GO J B5'(r), which 
a contradiction to assumption. So the proposition is true. 

Definition 4.7: (Necessary element set) NES(T)-{ oi : cxj is necessary element with respect to T}. 

Corollary 4.8: For any S € OES(T), NES(T) C E . 
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Algorithm 4 9- CNESI 




Tnnnt- d — (M V„\ M 




(Tj O Y i) h) T- the order 5] 




4 (Ti /To /T V * 

\ u xi u 2i * • • j u n J j 




Initialization: Set iV£5 := 0; 




AddElement: for i = 1 to ti do 




Uc^lll 11 Lx IS IIOL UlagllOSdUie 




with respect to £ \ {<Jj} and T 




then 




A/£S := NESU{aj}; 




end; 




Output: Return NES; 




Fig. 6. Algorithm:NES(T). 



Proo/.- For any S D G OES(T), and any cr, G NES(T), there exists <7j G S D (see Proposition l4~6l So 
NES(T) C S . ■ 
We introduce NES(T) to reduce the computing time. We partition the finite events space £ into two 
disjoint subsets NES(T) and S \ NES(T). The set NES(T) must include all the elements in OES(T). 
If we get NES(T) first, Algorithm @3]in this section need only compute in set £ \ NES(T). This may 
reduce computing complexity. 

D. Algorithm Complexity 

Suppose the time (of whether G is not diagnosable with respect to minOES and T) is Tq, where 

T G = 0{m 2 n). 

In Algorithm 12.11 in "Removal" recycle, the bad time is n x Tg, and the time-complexity of Algorithm 12.11 

is 0(m 2 n 2 ). 

In Algorithm |4.3| in "Initialization" recycle, it first calls Algorithm 12.11 to get minOES, the bad time 
is 0(m 2 n 2 ); and then it get set H, its a 0-1 pack problem, the bad time is 0(n x cmS); in "Testing 
diagnosability" recycle, the bad time is ng xTg; therefore the time-complexity of Algorithm 14.31 is 0(m 2 x 
n x ng). 

In Algorithm 14.91 in "AddElement" recycle, the bad time is nxTc, so the time-complexity of Algorithm 14.91 

is 0(m 2 n 2 ). 

Because S a C £, n = |S| in this section is greater than n = |S D | in Section IIII-BI 
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(71,04 01,03 




02,03 02,04 
Fig. 7. An Automata. 

V. Examples 

A. Example of Algorithm in Section \III-B\ 

Let us consider the system which is visualized as Fig. [7J 

From Fig. [71 we know Q = {90,91,92,93,94} and S = {01,02,03,04}- It is easy to compute that: 

<2(ci) = {qi,q2}, (3(02) = {93,94}, Q(<r 3 ) = {92,93}, Q(^i) = {90,91,94}- 

Diagnosability of the circuit depends on S and T. Let the desired partition T = {{90}, {91}, {93}}- We 
consider the following two examples for S G . 

Let S Q = {01,02}, we first use Algorithm 13.31 to compute A T a . In "Initialization" section, set m = 
5, n = 2, so = s\ = S2 = S3 = S4 = 0. In "Intersection"recycle, when step i = 1(01), we get si = S2 = 1; 
when step i = 2(02), we get S3 = S4 = 2. The final result is sq = 0, si = S2 = 1, S3 = S4 = 2. And then we 
send the result to Algorithm 13.41 By using quicksort algorithm, we get the result so < si = S2 < S3 = S4. 
In the last statements of Algorithm 13.41 we find si = S2- But in the desired partition T, q\ € {91}, and 
92 £={91}, so we get the OFD = FALSE in final. So the system is not diagnosable with respect to T and 
S . 

Let S D = {01,02,03}, we first use Algorithm 13.31 to compute /\ T a . In "Initialization" section, set 
m = 5, n = 3, so = si = S2 = S3 = S4 = 0. In "Intersection"recycle, when step i = 1(01), we get 
si = S2 = 1; when step i = 2(02), we get S3 = S4 = 2; when step i = 3(03), we get S2 = 5, S3 = 6. The 
final result is so = 0, si = 1, S2 = 5, S3 = 6, S4 = 2. And then we send the result to Algorithm 13 .41 By using 
quicksort algorithm, we get the result sq < si < S4 < S2 < S3. In the last statements of Algorithm 13.41 all 
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the values of Sj(j = 0, 1, 2, 3, 4) are not equal, and we get the OFD = TRUE in final. So the system is 
diagnosable with respect to T and E G . 

B. Lost of Algorithm 12.71 
Let 

M = {<Xi,cr 2 , . . . ,0io}. 

And the cost of M is 

C(M) = {27, 23, 20, 15, 10, 9, 7, 5, 4, 1}, 

G is diagnosable with respect to given T and an element of OES{T). 

OES{T) = {Ei = {a 3 ,a 5 , a 7 , a 10 }} U {E 2 = {a 3 , a 5 , a 8 , a 9 , a 10 }} U {£ C M : E x C S, or, £ 2 C S}. 
Using Algorithm [27Q we can get the minOES = E 2 , but the cost of E 2 is 20 + 10 + 5 + 4 + 1 = 40. The 
cost of Ei is 20 + 10 + 7 + 1 = 38, so the minOES is not the minimal cost of OES(T). 
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C. Example of Algorithm \4.3\ 

Using Algorithm 14.31 in "Initialization" section, we get minOES = {0-3, 0-5, as, 09, 010}, ImS = minL(min 
3, cmS = C{minOES) = 40. And then we get 

{(71, 0-5,(Tlo}, 
{0-l,Cr 6 , CT 9 }, 

{o"i,cr 6 ,aio}, 

{<Ti,o- 7 ,cr 8 ,<7io}, 

{0-1,0-7, Og, CT10}, 

{01, a 8 , C9,<7io}, 

{a2,a4,a w }, 

{a 2 ,a 5 ,a 7 }, 

{a 2 ,a 5 ,a 8 ,a w }, 

{o 2 ,o 5 ,o Q ,oi }, 

{a 2 ,aQ,a 7 ,ai }, 

{a"2,cr6,cr 8 ,CJio}, 

H = < 

{o-2,cr 6 ,cr9,CJio}, 
{a 2 , cry, cr 8 , 0-9,0-10}, 

{o- 3 ,o- 4 ,o- 8 }, 

{0-3, 0-4, O-g, O-10}, 
{0-3,0-5,0-6,0-10}, 
{0-3, 0-5, 0-7, O-10}, 
{0-3, 0-5, 0-8,0-9, O-10}, 
{0-3, 0- 6 ,0- 7 , 0-9}, 
{0-3, 0-6,0-7, O-10}, 
{0-3, 0- 6 ,0- 8 ,0- 9 , O-10}, 
{0-3, 0-7, 0-8,0-9, O-10}, 

other not empty subset of above set. 

In "Testing diagnosability" recycle, we find that only two elements of H(Ei and E 2 ) are diagnosable, 
and C(Ex) < C(E 2 ). Therefor we get the minimal cost of OES{T) is C(Ei). 

If we consider the set NES. From Algorithm 14.91 we get the set NES. We partition the finite events 
set £ into disjoint subsets NES = {0-3, 0-5, 010} and £ \ NES = {0-1,0-2, 0-4, 0-6,0-7, 0-8,0-9}. Now in Algo- 
rithm [43J we use the set (E\NES) as the set S. The computing procedure is as follows: in "Initialization" 
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H = < 



section, we get minOES = {erg, erg}, ImS = minL(minOES) = 8, cmS = C(minOES) = 9. And then 
we get that 

We}, 

{0-8,0-9}, 

other not empty subset of above set. 
In "Testing diagnosability" recycle, we find that only two elements of #({07} and {as, erg}) are diagnos- 
able, and C({a 7 }) < C{{a 8 , a 9 }). Hence we get that the minimal cost of OES{T) is C({a 7 }) + C(NES). 
The result is the same as that by using the method above, but the complexity is greatly reduced. 

VI. Conclusion 

In terms of some problems in off-line diagnostics [15], in this paper, we present some off-line diagnostic 
algorithms to overcome the shortcomings. We give a general method of judging a system's off-line diagnos- 
ability, which is a polynomial-time algorithm. And we give an algorithm of how to find the minimal set in 
all observable event sets. Of course, another issue worthy of further consideration is the on-line diagnostic 
algorithms of the minimal cost in DESs. We would like to consider it in subsequent work. 
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